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The logarithmic violations of the area law, i.e. an "area law" with logarithmic correction of the 
form S ~ L'^~^ log I/, for entanglement entropy are found in both ID gapless fermionic systems with 
Fermi points and for high dimensional free fermions. The purpose of this work is to show that both 
violations are of the same origin, and in the presence of Fermi liquid interactions such behavior 
persists for 2D fermion systems. In this paper we first consider the entanglement entropy of a toy 
model, namely a set of decoupled ID chains of free spinless fermions, to relate both violations in 
an intuitive way. We then use multi-dimensional bosonization to re-derive the formula by Gioev 
and Klich [Phys. Rev. Lett. 96, 100503 (2006)] for free fermions through a low-energy effective 
Hamiltonian, and explicitly show the logarithmic corrections to the area law in both cases share 
the same origin: the discontinuity at the Fermi surface (points). In the presence of Fermi liquid 
(forward scattering) interactions, the bosonized theory remains quadratic in terms of the original 
local degrees of freedom, and after regularizing the theory with a mass term we are able to calculate 
the entanglement entropy perturbatively up to second order in powers of the coupling parameter for 
a special geometry via the replica trick. We show that these interactions do not change the leading 
scaling behavior for the entanglement entropy of a Fermi liquid. At higher orders, we argue that 
this should remain true through a scaling analysis. 



I. INTRODUCTION 

The study of entanglement, which is one of the most 
fundamental aspects of quantum mechanics, has lead to 
and is still leading to much important progress and ap- 
plications in different fields of modern physics such as 
quantum information[l], condensed matter physics[2-4], 
etc.. To name a few, it has lead to better undertanding 
of density matrix renormalization group (DMRG)[ - ]; 
it has also been proposed to be a tool for the characteri- 
zation of certain topological phases[>— j ()]. 

Among various ways of quantifying entanglement, in 
condensed matter or many-body physics efforts have 
mainly focused on the bipartite block entanglement 
entropy (von Neumann entropy) and its generaliza- 
tions (Renyi or Tsallis entropy). It has become in- 
creasingly useful in characterizing phases[ri] and phase 
transitions[12, 13]. The area law[14] is one of the most 
important results on entanglement entropy: it states that 
the entanglement entropy is proportional to the area of 
the surface separating two subsystems. However, thus far 
there are two important classes of systems that violate 
the area law: in gapless one dimensional (ID) systems, a 
logarithmic divergence [13, 15] is found where according 
to the area law the entanglement entropy should saturate 
as the size of the subsystem grows; in higher dimensions, 
for free fermions the area law is found to be corrected by 
a similar logarithmic factor log L[l(i-2.')], where L is the 
linear dimension of the subsystem. 

In this work, we first show that the scaling behavior 
of the entanglement entropy for systems with a Fermi 



surface is the same as that of ID systems with Fermi 
points[24-27]. We then seek for a generalization of the 
latter to interacting fermions in the Fermi liquid phase. 
We first develop an intuitive understanding via a toy 
model, showing that in this model the entanglement en- 
tropy has the same form as that given in by Gioev and 
Klich (GK) in Ref. [ ]. We then develop a more gen- 
eral and formal treatment using the method of high- 
dimensional bosonization[28-32]. This approach will not 
only lead to a reproduction of the result for free fermions 
obtained by GK based on Widom's conjecture['^ ', ? !], 
but will also lend itself to the inclusion and subsequent 
treatment of Fermi liquid type (forward scattering) in- 
teractions. 



This paper is organized as follows. In Sec. II, we de- 
scribe the toy model for which the entanglement entropy 
can be written in the same form as the GK result. Then 
in Sec. Ill, we briefly introduce the tool box of multi- 
dimensional bosonization, and apply it to free fermions to 
reproduce the GK formula. The main results of this work 
are presented in Sec. IV in which we calculate the entan- 
glement entropy of a Fermi liquid for a special geometry 
using a combination of multi-dimensional bosonization 
and the replica trick. We subsequently summarize and 
discuss our results. Some technical details are discussed 
in two appendices. 
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FIG. 1: (Color online) The toy model both in real space and 
in momentum space, (a) A set of parallel decoupled Id chains 
of spinless free fermions (dash lines); the subsystem division 
is represented by the solid lines, both convex and concave 
geometries, (b) Fermi surfaces of the toy model. 



II. THE INTUITIVE PICTURE - A TOY MODEL 

Consider a set of decoupled parallel ID chains of non- 
interacting spinless fermions with spacing a as shown in 
Fig. 1(a). Here we only consider d = 2 for simplicity, 
but this toy model is viable in general d dimensions. The 
asymptotic behavior of entanglement entropy in large 
L limit of a convex subsystem A of this model can be 
obtained by simply counting the number of chains that 
intersect A, and each segment contributes a (1/3) log L 
where L is the linear dimension of the subsystem[14, 35]. 
Due to the logarithm, different shapes only lead to dif- 
ferences at the area law level. Since each segment must 
have two intersections, we can count the intersections in- 
stead, which also automatically takes care of non-convex 
geometries. Although there is an additional correction 
for multiple intervals on a single chain[37-39], as long as 
only the log L behavior is concerned, that contribution 
is negligible. For L large enough, we can write the num- 
ber of these intersections as an integral over the surface 
of A projected onto the direction perpendicular to the 
chains times one half of the chain density, 1/a. To make 
contact with the GK result, we note this model also has 
Fermi surfaces as shown in Fig. 1(b) with a total "area" 



of 47r/a. This enables us to replace the density of chains 
by an integral over the Fermi surfaces of the system 
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where F indicates the occupied area in momentum space 
so its boundary dT is the Fermi surface(s). Therefore we 
can write the entanglement entropy as 
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where h is the direction along the chains which is also 
normal to the Fermi surface, and an overall factor of | 
accounts for the double counting of chain segments. In 
Eq. (1) we recover the GK formula in this special case 
but written in a slightly different way. In Ref. [](>] the 
entanglement entropy is given as: 
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where the real space surface integral is carried out over 
the subsystem whose volume is normalized to 1. The 
surface area is factored out as L'^~^ . However, in our 
formula the surface area ~ L'^~^ is implicitly included in 
the integral over the surface of the subsystem. 

We note that the model discussed in Ref.[''^'] is equiv- 
alent with our toy model, but motivated from a different 
perspective. In Ref.[2(i], models are constructed from the 
momentum space, either with Fermi surfaces as our toy 
model, or a square Fermi surface, and a boxlike and a 
spherical geometry are discussed. In contrast, our toy 
model is constructed from a real space perspective, and 
general single connected geometries are discussed. 

Motivated by the toy model, in this work we extend 
this intuitive understanding of GK's result to generic free 
Fermi systems and generalize it to include Fermi liquid 
interactions in two dimensions (2D) via high dimensional 
bosonization. Using the method of multi-dimensional 
bosonization, the Fermi liquid theory can be written as a 
tensor product of low-energy effective theories of quasi- 
ID systems similar to this toy model, along all directions. 
This provides us with a tool to treat the entanglement 
entropy of fermions in high dimensions, even in the pres- 
ence of interactions. 

At this point, we could also include forward scattering 
for each chain, and from ID bosonization we know that 
for spinless fermions this only leads to renormalization of 
the Fermi velocity, thus does not change the logarithmic 
scaling of the entanglement entropy for this toy model. 
This hints that the same conclusion might hold for Fermi 
liquids, as we can include Fermi liquid interactions in a 
similar way via high dimensional bosonization. Although 
as we show later, this is indeed true at the leading order, 
the situation is more delicate than it seems to be. The 
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Fermi liquid interactions couple a family of "toy models" 
aligned along different directions in the language of high 
dimensional bosonization, and lead to a correction to the 
entanglement entropy ^ 0{1) x \ogL. 




III. MULTI-DIMENSIONAL BOSONIZATION 



The scheme of multi-dimensional bosonization was first 
introduced by Haldane[2(S], followed by others[29-32]. 
The basic idea is to start with a low energy effective 
Hamiltonian (obtained through a renormalization group 
(RG) approach) restricted to within a thin shell of thick- 
ness A around the Fermi surface, kp — X/2 < \k\ < 
kF + )^/2. Then one divides this thin shell into N patches 
with dimensionality ~ A''"^ x A as shown in Fig. (2) in 
such a way that A ^ A ^ fcj? and A'^/kp «C A, where 
d = 2, 3 is the space dimension, A is the linear dimen- 
sion of the tangential extent of each patch. The condi- 
tion A ^ A minimizes inter-patch scattering; A <ti kp 
and A^ /kp <C A together makes the curvature of the 
Fermi surface negligible. In the end we shall take the 
limit A/kp — >■ 0, so that the sum over all patches can be 
converted to an integral over the Fermi surface. In this 
work, we treat the free theory in general d dimensions, 
but shall restrict ourselves to d = 2 when interactions are 
included. For an arbitrary patch labeled by the Fermi 
momentum ks at the center of the patch, we introduce 
the patch fermion field operator 




(a) 




i!{S- x) = e''=«-=^ 0{S; p)e'(P^'='^)-^ Vp, (3) 
p 

where tpp is the usual fermion field in momentum space, 

^ 1 1 if p lies in the patch S, 

' I if p lies outside patch S. 

The effective Fermi liquid Hamiltonian can be written as 

+ I d'^xd'^vYv{S,T-x-y)xl;\S]x)il}{S;x) W 

y.^\T-y)^{T-y), 

with m* being the effective mass, V{S, T; x—y) the effec- 
tive interaction in the forward scattering channels. Even 
though this model is restricted to special interactions of 
this form, forward scattering is known to be the only 
marginal interaction in RG analysis [40]. As the leading 
order contribution of the entanglement entropy is dom- 
inated by the low energy modes around the Fermi sur- 
face, it is sufficient to consider this model. Similar to 
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FIG. 2: Patching of the Fermi surface. The low energy theory 
is restricted to within a thin shell about the Fermi surface 
with a thickness A <C kp, in the sense of renormalization. 
The thin shell is further divided into A'^ different patches; 
each has a transverse dimension A"*"^ where d — 2, 3 is the 
space dimensions. The dimensions of the patch satisfy three 
conditions: (1) A ^ A minimizes inter-patch scattering; (2) 
A ^ fcf and A^/kp <C A together makes the curvature of 
the Fermi surface negligible, (a): Division of a 2D Fermi 
surface into A'^ patches. Patch S is characterized by the Fermi 
momentum ks- (b): A patch for d = 3. The patch has a 
thickness A along the normal direction and a width A along 
the transverse direction (s). 

the ID case, the bosonic degrees of freedom are the den- 
sity modes of the system, in this case defined within each 
patch of the Fermi surface: 

k 

(5) 

Though q is not explicitly bounded in the above defini- 
tion of the patch density operator, its transverse compo- 
nents qs± = (gsl, • ■ • , Is'h ■ ■ ■ ' 9s± (those parallel to 
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the Fermi surface) are limited q^^l^ G (—A, A) due to the 
patch confinement. Their commutation relation is 



where 



n 



[J{S; q), J{T; p)] ^ Ss,tS'^+p,o E ^(■^^ 

k 

X [e{S-k-q)-e{S;k + q)]nk 
= <5s,T<5^+p,of^ (^s ■ (D02{qsA.) + 0(A/A), 



(6) 
(7) 

(8) 



A [Lo/(27r)] , Uk — (V'fcV-'fc) is the occupation 
number of state with momentum fe, hg is the outward 
normal direction of patch S, qs± represents all other 
component(s) of q that are perpendicular to ng, and Lq 
is the linear dimension of the entire system. The appear- 
ance of (5q+p,o is a result of momentum conservation. The 
calculation of the commutator is reduced to computing 
the difference of occupied states, i.e. the area difference 
below the Fermi surface, between the two 9 functions 
ie{S; k-q)- e{S; k-p)) as indicated by Eq. (6). This 
is similar to ID bosonization. If we consider both k and 
q to be ID momenta, Eq. (6) would give us the ID 
bosonization commutator. The 2D result Eq. (7) is sim- 
ilar, because the Fermi surface confined within the patch 
is essentially flat thus the dispersion is ID. That leads to 
the fis ■ q dependence of the commutator as that of the 
ID case, even for qg^^ ^ 0. The difference is that, as il- 
lustrated in Fig. (3), due to the patch confinement on the 
transverse direction(s), when qsj. increases k ± q would 
increasingly find itself outside the patch thus not con- 
tributing to the commutator. According to Fig. (3), one 
can see that this gives rise to the factor 02{qs±), which 
diminishes the commutator at large qs± ■ It is usually ne- 
glected in literature because the long wavelength limit is 
taken[29-32]. However, as it is important in the present 
context to correctly count the number of total degrees of 
freedom, this d2{qs±) factor cannot be neglected because 
it comes from counting the transverse degrees of freedom. 
To simplify things, we replace 02(qs_l) by 



e2iqs±) = 1 for -A/2 < < A/2, 



(9) 



and we also limit qs± to this range. This approximation 
makes it easier to do Fourier transform while keeping the 
total degrees of freedom intact. To see that, it is suffi- 
cient to consider one direction, comparing the area en- 
closed by the two different functions; 92{q±) = 1 — q±/A 
over the range (—A, A) and 6*2 = 1 over the range 
(—A/2, A/2). Both functions enclose the same area thus 
the same number of states. This approximation can also 
be interpreted as relaxation of the hard wall cutoff in Eq. 
(5), softening of the step function 9{S; k). In Eq. (5), q 
is not bounded while k is bounded by 9{S; k). If we relax 
the restriction on k on the transverse direction, allowing 
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FIG. 3: (Color online) Origin of the bosonic commutator of 
patch density operators illustrated for d — 2. As shown in Eq. 
(7), the commutator is reduced to computing the difference 
of occupied states , i.e. the area difference below the Fermi 
surface, between the two 6 functions {6{S; k~q)-'9{S; k+p)). 
The solid box indicates the original patch, or 9{S; k). The red 
line shows the Fermi surface. Both 6[S; k — q) and 6[S; k~\-p) 
are denoted by dashed boxes. The occupied part in 6{S\ k+q) 
is denoted by blue, that of 6{S; k — q) is denoted by red, and 
the overlapping region is denoted by yellow. Subtracting the 
remaining blue area from the red, we obtain that 6{S;k — 
q) occupies (A — qs±)qns more states, which gives us the 
commutator. 



k with > A/2 in the summation, but require qgj_ 

to be bounded within the patch, we would obtain the 
alternative 92{q±). 

Using from now on the above approximation, we con- 
struct the local bosonic degrees of freedom (t>{S; x) = 
(j){S;xs,xs±) as 



J{S;x) ^ Vfid.j;s(l){S;xs,xs±), 



(10) 



where J{S;x) = '^^e^'^'^ J{S\q), xs — x ■ fis, and 
xs± — X — [x ■ ns)ns- The commutation relations for 
the (b's are then 



[d.,,<f>iS; x),^{T; y)] - i27rn6s.TS{xs - ys) 



n 



'sin(A(4"|-y^")))^ 



(11) 



ysl) 



which is the bosonic commutation relation we are look- 
ing for. The factor JJa i ) ^I'ising from 

transverse directions must be treated with care in differ- 
ent circumstances. In most literature, the focus is the 
physics at large length scale I ^ 1/A; therefore, this fac- 
tor is usually approximated by 8'^^^(xs\_ — ys±) which 
is good in that limit without further discussion. This is 
also what we shall do for most of the time unless noted 
otherwise: 



i2TTnSsJ,S{xs-ys)S'' '(xs^-ys^ 



d-l/ 



(12) 
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However, the more accurate expression (11) is useful for 
us to understand how to count the transverse degrees of 
freedom correctly. It tells us that the transverse degrees 
are not independent on the short length scale I < 1/A. 
More importantly, later on we need to consider the limit 
(5'*"^(a;s_L - ys±)\ys±^xs±'- without Eq. (11), this limit 
would be ill-defined. 

With the above, the Hamiltonian H[ip^,ijj] is found to 
be quadratic in terms of these J{S; q)'s: 



JiS;-q)JiT;q) 



S,T;q -- (13) 

+V{S,T;q)J{S;-q)J{T;q), 

where V{S, T; q) is the Fourier transform of V'(iS', T\x — 
y) . So it is also quadratic in the bosonic fields associated 
with the J{S;qys. 

A. Entanglement Entropy of Free Fermions 

The kinetic energy part of Eq. (4) or its bosonized 
version Eq. (13) can be written in terms of the boson 
fields constructed above as: 



S;q 



2t:vI 



nv 



J2 [ d'x{d,,<j>{S;x)f 



(14) 



We see that there is no coupling between different 
patches. The theory is thus formally a tensor product 
of many independent theories, one for each patch. We 
can therefore calculate the entanglement entropy patch 
by patch and sum up contributions from each patch in 
the end. Within a single patch there is no dynamics in 
the perpendicular direction as dictated by the Hamilto- 
nian, and the problem is reduced to a one dimensional 
problem! Note that transverse degrees of freedom are 
not completely independent. According to Eq. (11), the 
commutator is non-vanishing for xs± 7^ ys± up to a 
length scale ~ 27r/A. This is a consequence of restricting 
qs± to within the range [—A/ 2, A/ 2]. Physically one can 
view this as discretization along the transverse direction 
due to a restricted momentum range, similar to the re- 
lation between a lattice and its Brillouin zone. In this 
view the single patch problem is reduced to a ID prob- 
lem with a chain density of (A/(27r))'* ^. Therefore, the 
Hamiltonian (14) becomes 



Hn = 



nv 



E 



dxs {dxs4>{S;x)Y 



(15) 



Note that the bosonized theory of a single patch is chi- 
ral. To directly make use of our toy model, we need to 
consider two patches having opposite hs simultaneously. 



This is because for a ID fermion model at non-zero filling, 
there are two Fermi points. Both need to be considered 
to construct well-defined local degrees of freedom. Once 
we consider such two patches together, it is more conve- 
nient to combine the two chiral theories into a non-chiral 
theories. This is also what we will do for the rest of this 
work. Introduce the non-chiral fields 



x{S;x) 



V2 



-S;x)), 
-S;x)), 



(16) 



where —S indicates the patch with normal direction op- 
posite to that of patch S: fi-s = —f>-s- One finds that 
X and if are mutually dual fields with S restricted to one 
hemisphere, but dxs^^ and tp now commute while x £^nd 
if have a non-trivial commutator: 

MS; x),dyMS; y)] = [x{S; x),dy,xiS; y)] = 0, 

[d^MS; x),x{T; y)] = [d.,sX{S; x), <f{T; y)] (17) 
= 2iTTnSs.,TS{xs - ys)S''^^{xs± ~ ys±)- 

Therefore, two patches with opposite fis are equivalent 
to a set of ordinary ID boson fields. Throughout the 
rest of this work, we shall assume this chiral-to-nonchiral 
transformation is done, and when we refer to patches 
we always refer to the two companion patches that form 
a non-chiral patch together. For the non-chiral boson 
theory, it is known that the entanglement entropy of a 
single interval (with two end points) is (1/3) log L. 

Before we proceed further, we note that the relation 
between boson fields and the original fermion fields is 
not completely local. However, the underlying physical 
quantity that matters is not the fields, but the fermion 
density, or in other words, the fermion number basis one 
chooses to expand the Hilbert space of the problem. This 
physical basis is also what one uses to do the partial trace. 
It is known that the fermion density operator obeys a 
locally one-to-one corresponding relation to the boson 
fields. Thus we argue that in ID the nonlocal relation 
between the fermion and boson fields does not affect the 
partial trace operation, so as the calculation of entangle- 
ment entropy. 

By referring to our result for the toy model, the con- 
tribution from a single patch is readily given 



S{S) ^^\og L£jns-dS,\x (^A^ 



(18) 



where an additional factor of 1/2 has been introduced 
in order to count only once each pair of patches forming 
a non-chiral theory. Identifying nsA'^^^ as the surface 
element at the Fermi surface dSk and taking the N 00 
limit, the total entanglement entropy is 

S = log^ / / \dSk ■ dS^l (19) 

12(27r)'' JdAJdT 

So we recover the GK result for generic free fermions. 
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B. Solution for the Fermi liquid case and 
non-locality of the Bogoliubov fields 

When Fermi liquid interactions (forward scattering) 
are included, the full Hamiltonian will no longer be di- 
agonal in the patch index S. But it is still quadratic 
in terms of the patch density operators, i.e. the bosonic 
degrees of freedom, and can be diagonalized by a Bo- 
goliubov transformation. According to Eq. (7) and ig- 
noring terms of ©(A/A), one can define a set of boson 
creation/annihilation operators a^{q)/a{q) as follows: 



q,nsq>Q 



v4 



ns ■ q\ 



It can be shown that the full Hamiltonian is diago- 
nal in q, and it can be diagonalized by a Bogoliubov 
transformation[.)2] independently for each q sector. In 
Ref. [:!■-], only a Hubbard- C/ like interaction is considered 
for practical reasons. But in principle, such a Bogoliubov 
transformation also applies to general interactions: 



J 



(21) 



where f{S, l;x — y) is typically long-range, even for the 
short-range Hubbard-J7 interaction. For more general 
cases, with further q-dependence in the interaction, the 
non-locality would only be enhanced. The loss of locality 
prevents us from calculating the entanglement entropy 
directly using those eigen modes, since it is difficult to im- 
plement the partial trace using those non-local degrees of 
freedoms. Therefore, although the Bogoliubov fields have 
a local core as we would expect for Fermi liquids from adi- 
abaticity, they do acquire a nonlocal dressing due to inter- 
action. Though in principle the partial trace can be done 
with those Bogoliubov fields, such nonlocality makes it 
difficult and we have not been able to do it, which further 
renders calculating the entanglement entropy impossible. 
This is very different from the ID theory, where for lo- 
cal interactions the eigen fields remain local, since there 
are only two Fermi points. There the transformation can 
never involve such angular q-dependence due to limited 
dimensionality. Despite these technical difficulties, the 
non-locality may suggest possible corrections to the en- 
tanglement entropy. This is indeed the case as revealed 
by our later calculation for Fermi liquid interactions, al- 
though in this case such extra contributions are only of 
0(1) X logL which is of ©(l/L) comparing to the leading 
term. This shows that the mode-counting argument in 
Ref. ] , though correctly suggesting the log L violation 
to the area law for Fermi liquids, does not always fully 
account for all sources of entanglement entropy. 



where both i and j refer to the patch index, aj and /3j 
are the Bogoliubov bosonic annihilation operators that 
diagonalize the Hamiltonian. With proper choice of u's 
and z;'s, the Hamiltonian is readily diagonalized. Ref. 
[32] solves the Hubbard-C/ like interaction and provides a 
successful description of Fermi liquids, even in the strong 
U limit. 

However, even for this simple case in which U has no q- 
dependence, the Bogoliubov transformation still depends 
on q. To be more precise, Uij and Vij will depend only 
on the angle between the patch normal direction ris and 
q, leading to discontinuities in the derivatives at q = 
0. Consequently, the real space fields constructed from 
the Bogoliubov operators aj and Pj are no longer local 
with respect to the original boson fields. The real space 
Bogoliubov fields are constructed in a manner similar to 
Eq. (20): 



q,nsq>Q 



a^{S;q)e-'''-'' - a{S;q)e"i '- 
Vns ■ q 



(22) 

Then one can show that the original local degrees of free- 
dom 0(5'; x) can be expressed in terms of above Bogoli- 
ubov fields as 

^{S; x) = 4>{S; x)+ I dy^ fi^, h^- y)0(/; y), 

I 

(23) 



IV. ENTANGLEMENT ENTROPY FROM THE 
GREEN'S FUNCTION 

In order to preserve locality, we need to work with the 
original local degrees of freedom. To do that, we adopt 
the approach used by Calabrese and Cardy[41] (CC) on 
calculating the entanglement entropy of a free massive 
ID bosonic field theory. The calculation is done in terms 
of the Green's function by applying the replica trick. In 
our case, we find that the CC approach can be generalized 
in a special geometry for solving the interacting theory 
which is quadratic after bosonization. In this way, we 
avoid diagonalizing the Hamiltonian and thus the nonlo- 
cality issue. However, we do have to regularize the theory 
by adding a mass term by hand. In the end we shall take 
the small mass limit, and replace the divergent corre- 
lation length ^ ^ I/to by the subsystem size L. The 
regularization procedure facilitates the calculation, but 
also strictly restricts us to computing the entanglement 
entropy only at the logi level. 

In this section, by using the replica trick we convey 
the calculation of entanglement entropy into computing 
the Green's function on an n-sheeted replica manifold. 
We first demonstrate the method by applying it to free 
fermion theory in d-dimensional space; then based on it 
we compute the entanglement entropy perturbatively for 
a simple Fermi liquid theory in powers of the interaction 
strength up to the second order. 



A. The Replica Trick and Application to ID Free 
Bosonic Theory 

In this part, we briefly describe the rephca trick in 
(1 + 1) space-time dimensions {{l + l)d) so that later on 
we can straightforwardly generalize it to (2 + 1) space- 
time dimensions ((2 -I- l)d) accordingly for our problem. 

The replica trick makes use of the following identity: 

Sa = -tr (pa In pa) = - lim ^^tr p^. (24) 

n-i-i on 

To compute tr p^, CC use path integral to express the 
density matrix p in terms of the boson fields 

p({0(x)}|{0(x')'}) = Z-i({0(x)}|e-^|{0(x')'}), (25) 

where Z = tr e~^^ is the partition function, (5 is the 
inverse temperature, and {(p^x)} are the corresponding 
eigenstates of (t){x): (j){x)\{(j){x')}) = (j){x')\{(j){x')}) . p 
can be expressed as a (Euclidean) path integral: 



p=Z-^ I mx, t)] n S{cl){x, 0) - 0(x)') 
X J|5((/.(x,/3)-0(a;)")e-^^ 



(26) 



where Se = LeCLt, with Le being the Euclidean La- 
grangian. The normalization factor Z, i.e. the partition 
function is found by setting {cf){x)"} — {4>{x)'} and inte- 
grating over these variables. This has the effect of sewing 
together the edges along t — Q and t — (5 to form a 
cylinder of circumference (3 as illustrated in Fig. (4) (left 
panel) . 

The reduced density matrix of an interval A = {xi, x f) 
can be obtained by sewing together only those points 
which are not in the interval A. This has the effect of 
leaving an open cut along the line r = which is shown 
in Fig. 4 (right panel). To compute p^, we make n copies 
of above set-up labeled by an integer k with 1 < k < n, 
and sew them together cyclically along the open cut so 
that 0(x)'j, = 0(x)';!^Jand (/)(a;)'„ = 0(a;)"] for aU x e A. 
In Fig. 5(a) we show the case n = 2. Let us denote 
the path integral on this n-sheeted structure (known as 
n-sheeted Riemann surface) by Zn{A). Then 



tr pI 



so that 



Sa 



ZniA) 



,. d Z„iA) 
lim . 

n^i dn Z" 



(27) 



(28) 







( 






FIG. 4: (Color online) Path integral representation of the 
reduced density matrix. Left: When we sew (j}{xY = cj>{x)" 
together for all x's, we get the partition function Z. Right: 
When only sew x ^ A together, we get pA- 




(a) 




(b) 



FIG. 5: (Color online) Formation of the n— sheeted Riemann 
surface in the replica trick. By sewing n copies of the re- 
duced density matrices together, one obtains the replica par- 
tition function Z„. In the zero temperature limit, /3 — ^ oo, 
each cylinder representing one copy of pA becomes an infinite 
plane. Those n— planes sewed together form a n— sheeted Rie- 
mann surface in Fig. 5(b) which can be simply realized by 
enforcing a 2n7r periodicity on the angular variable of the po- 
lar coordinates of the (1 -I- l)d plane instead of the usual 27r 
one. (a): n copies of the reduced density matrices. For clar- 
ity only n = 2 is shown, (b): Visualization of a n— sheeted 
Riemann surface. 



If we consider the theory as that of one field living 
on this complex n— sheeted Riemann surface instead of 
a theory of n copies, it is possible to remove the replica 
index n from the fields, and instead consider a problem 
defined on such an n-sheeted Riemann surface which can 
be realized by imposing proper boundary conditions. 

In Ref.[4i], CC consider the entanglement entropy be- 
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tween the two semi-infinite ID system (i.e. cutting an 
infinite chain into two halves at a; = 0) for free mas- 
sive boson fields. For such geometry, as illustrated in 
Fig. 5(b), the n— sheeted Riemann surface constraint is 
realized by imposing a 2n7r periodicity on the angular 
variable of the polar coordinates of the (1 -I- l)d plane 
instead of the usual 2tt one. In this way, the (1 + l)d 
variable x = (x, r) acquires n branches a;„, and each 
branch corresponds to one copy of (j). Notation-wise this 
corresponds to 

<j){x, T)k (pixk) (t)[x), (29) 

and the sewing conditions 0(a;)'^. ~ '/'(a;)fc+i simply be- 
comes the continuity condition for 4'{x) across its con- 
secutive branches. Here we use a generalized polar coor- 
dinate: X = (r, 9) with < r < oo, and < 6* < 2rm. 

The massive free boson theory considered by CC is 
defined by the following action 

The (1 -I- \)d bosonic Green's function G^^l{r,r') — 
{(j){r)4>{r')) on the n— sheeted Riemann surface satisfies 
the differential equation 

(-V^ + m2)Gg=^(r-r'). 

To compute the partition function, one can make use of 
the identity 

^logZ„ = - d^+'xG^-\x,x). (30) 

Note that here the integration is over the entire 
n— sheeted space. The above is applicable to general 
quadratic theories of bosons, and will be applied by us 
later to bosonized theories of interacting fermions. Here 
we use G'^"^(a;, a;'), a general two point correlation func- 
tion on the rt— sheeted Riemann surface in d-dimcnsional 
space for later use, instead of the specific Gq"^'' defined 
above. Accordingly, Sa is then given as 

Sa^- lim Ae-5/'i™'/'*°+'^(G<"'(a3,x)-nG(i'(x,a,))^ 

n->-i dn 

(31) 

Here and in the following, we will leave it understood that 
the first term in the integrand is integrated over the n- 
sheeted geometry, whereas the second is integrated over 
a one-sheeted geometry. There should be no confusion as 
the superscript of G generally indicates the geometry. 

The benefit of the above approach is that the two point 
correlation function or Green's function, defined in terms 
of certain differential equation obtained from the equa- 
tion of motion, can be solved for on the n— sheeted Rie- 
mann surface thus enabling us to compute the entangle- 
ment entropy. Although CC's work only considers mas- 
sive (1 -I- \)d boson fields, it is also applicable to our 



case. The price one has to pay is to to introduce a mass 
term for regularization. At the end of the calculation 
the inverse mass, which is the correlation length of the 
system, shall be considered to be on the same scale as 
L: 1/m ^ L, where L is the characteristic length scale 
of the subsystem. The validity of such consideration is 
well-established in other cases, [3">, 42] where the corre- 
lation length is either set by finite temperature or mass. 
The only modification necessary to apply the above to a 
bosonized Fermi surface in higher dimensions is to intro- 
duce a sum over the patch index. 

B. Geometry and Replica Boundary Conditions 

Through the remainder of this work, instead of the gen- 
eral geometry considered before, we work with a special 
half-cylinder geometry as shown in Fig. 6(a): the sys- 
tem is infinite in the x direction while obeying periodic 
boundary condition along the y direction with length L. 
The system is cut along the y axis so that we are com- 
puting the entanglement entropy between the two half 
planes. We require L to be large so that it can be con- 
sidered ^ oo unless otherwise noted. 

We choose such this simple geometry for the following 
reasons. Cutting the system straight along the y direc- 
tion, yielding a two half-plane geometry, is a straight- 
forward (2 -|- l)d generalization of the semi-infinite chain 
geometry considered in CC. It makes any straight line 
intersect the boundary only once, dividing it into two 
semi-infinite segments, for all patch directions as in the 
ID case, except for lines parallel to the hs — y patch di- 
rection. The degrees of freedom associated with this spe- 
cial patch do not contribute to the entanglement entropy, 
since they are not coupled (have no dynamics) along x, 
and are of measure zero in the large patch number limit 
anyway. 

For this simple geometry, the (2 + l)d n-sheeted geom- 
etry is constructed from n identical copies 

S-" = {(a;,y,T) e M X M X M}, (32) 

sliced along "branch cuts" 

C" = {(a;,y,T) e X M X {0}}, (33) 

and then appropriately glued together along these cuts. 
This happens exactly as in ID, and the y coordinate is 
so far a mere spectator. This defines an n— sheeted, or 
in this case more appropriately the n— layered, replica 
manifold which is a simple enough generalization of the 
(1 -t- l)d case. The n-sheeted Riemann surface, as dis- 
cussed in Sec. IV A and shown in Fig. 5(b), now acquires 
an extra direction y perpendicular to the x — f plane. It 
can still be implemented by imposing the same 2n7r pe- 
riodicity boundary conditions on 9, the angular variable 
of the polar coordinates {x,t) = {rcos9,rsm9) in the 
X — f plane. Therefore, we can safely make use of the 
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sheet via 



(a) 




(b) 



FIG. 6: (Color online) The half-cylinder geometry and equiv- 
alence of boundary conditions in a; — f and ns — f planes. 
The system is infinite in the x direction while obeys peri- 
odic boundary condition along the y direction with length L. 
The system is cut along the y axis so that we are computing 
the entanglement entropy between the two half planes, (a); 
The half-cylinder geometry, (b): The projection of rs onto 
the X — T plane. Consider polar coordinates of an arbitrary 
fis — T' plane (the blue plane). Since the polar coordinates in 
the x — T plane satisfies the 2mi periodic boundary condition, 
consider the one-to-one projection of the vector rs onto the 
x — T plane. Consider, if we move the vector in the x — t 
plane around the origin n times (the red circle). Due to the 
one-to one mapping, rs should also move around the origin 
n times (the blue "circle", it is actually a eclipse), thus obeys 
the 2n7r periodicity as well. 



CC result, i.e. the solution to the Green's function on a 
n— sheeted Riemann surface, to the free fermion theory, 
and can further use it as a starting point for treating the 
interacting theory. This is obviously true for the patch 
with hs — X, but it also holds for general fis as we shall 
validate as the following. 

For a general patch direction ns, the noninteracting 
Green's function associated with this patch embodies cor- 
relations in the affine hg — f "planes". The geometry of 
each such "plane" is that of the n-sheeted Riemann sur- 
face of the (1 -I- l)d problem, as we will now argue. With 
each patch direction we thus associate a different foli- 
ation of the n-sheeted (2 + l)d geometry into (1 + l)d 
counterparts. 

To be more precise, for given patch S, instead of 
the Gartesian coordinates {x,y,T), we consider a par- 
allel/perpendicular decomposition {xs,xs_l,t) for each 



(x, y) = xsns + xs±ns± 



(34) 



where the hsi_ are the perpendicular unit vectors aligned 
with the patch S. The natural choice of coordinates for 
a given patch is to choose polar coordinates within the 
ns — T' plane: 



rs = (xs + xs± 



"-s_L 
n% 



t) ^ (rs cos 9s, rssin0s) , (35) 



because these are the coordinates in which the fis — f 
planes restricted to each sheet are naturally glued to- 
gether by extending the range of 0s to 27rrt, as we will 
now show. The shift xs±ngj^/hg of xs is necessary as 
to ensure that a; = 0, the location of the onset of the 
branch cut, corresponds to rg = which is what makes 
these coordinates convenient. The ns — f planes are now 
defined by fixed xs±- 

If we can establish that the 27rn periodicity of 9 is 
equivalent to a 27rn periodicity of 9s, then GG's solution 
would be justified in the above set-up so that the non- 
interacting Green's function Go{S, S; ,rs,9s,r' s,9' s) 
can be expressed through GG's result. This can be 
achieved by establishing a one-to-one correspondence 
(mapping) between 9 and 9s- The mapping is intuitively 
constructed, as shown in Fig. 6(b), as the vertical projec- 
tion from the fis — t plane (the blue plane in Fig. 6(b)) 
onto the x — t plane along the y direction. Gonsider mov- 
ing the projection of rs in the x — t plane around the 
origin n times (the red circle). It is clear that rs (on 
the blue ellipse) follows its projection while also mov- 
ing around the origin n times, always being on the same 
sheet. In particular, the branch cut is always traversed 
simultaneously for 9 = 9s = mod 27r. The fis — f 
planes, the leaves of our foliation, thus have the familiar 
l-|-lc? n-sheeted geometry, and 9s obeys the same 2n7r 
periodicity as 9. 

Finally, the periodicity condition of the y direction is 
necessary for the total entanglement entropy to be finite; 
it also provides the only length scale for the subsystem 
which is needed for extracting the scaling behavior of 
entanglement entropy. However, if we are only concerned 
with the integral form of the entanglement entropy as in 
Eq. (19), not requiring it to be finite as a whole, but 
rather requiring only the entanglement entropy per unit 
length to be finite, we may take the y direction to be 
infinite. This point of view will be taken here and in the 
following in order to simplify our calculation. 



C. Entanglement entropy of free fermions revisited 

In order to treat the interacting theory, in this section 
we re-derive the free fermion result for the half cylinder 
geometry via the replica trick. Later we shall general- 
ize the method to include interactions. Rewriting the 
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Hamiltonian Eq. (14) in terms of the non-chiral fields, 



and adding the mass term by hand, we have 



(36) 



r 



For convenience, we use the Lagrangian formahsm and 
work with the ^(S\ x) representation through the rest of 
this work. 

Switching to imaginary time t — )• ir, and rescahng the 
coordinates in the fohowing manner: 



V 



167r^t;|,fi m' 



/ ?LY X 



(37) 



we obtain the following Lagrangian density in the (p rep- 
resentation : 

9 

?77 

^ = -^[idMS;x)r + {ds^{S;x)f + MS-x)n 

(38) 

Then we can work out the Euler-Lagrangian (E-L) equa- 
tion of motion. Making use of the E-L equation of mo- 
tion, we find the Green's function Gq^\s,T;x,x') — 
{T(p{S;x)ip{T;y))Q satisfies the following differential 
equation: 



(92 + 92^-1) G^"^(5,T;a;,y) 

C6s,tS{t - Ty)S{xs - ys)S'^^^{xs± - ys±) 



(39) 



d-2 



where C = 27rf7m'*-i (y • The rescaling 

makes the Green's function dimensionless, thus easier to 
handle when it comes to computing J d'''^^xG^^^\x,x). 
The extra factor C generated on the right hand side (rhs) 
will be canceled by the Jacobian of the integral over the 
Green's function, leaving only a factor of l/m? . All that 
needs to be computed is then an integral over the di- 
mensionless G. Therefore, it is legitimate to ignore this 
factor from now on. The (5-functions originate from the 
commutator Eq. (11), and are coarse-grained. After we 
include the patch index, perform the integral over m^, 
and take the n derivative, Eq. (31) becomes 



Sa^\ logim'al) hm ^ V(Cg(S; n) - nCaiS; 1)), 
2 n-i-i on ^ 



where Oq is an ultraviolet cutoff, and 



CG{S\n)^ J d'^+^xG''^\s,S;x,x). 



(40) 



(41) 



The exponential factor in Eq. (31) becomes one after the 
n — > 1 limit is applied. Note that ^ log{m^a^) ^ — logL. 
Our major task is now computing CaiSin). 

Observing that there is no xs± dependence on the left 
hand side (Ihs) of Eq. (39), we can write 

G[,")(S, T; X, y) = Ss^tS^-^xsi^ - ys±)Gg(5; r^,ry), 

and we obtain a (1 + l)d equation 

-{dl+dl^ + l)G^^l{S-rs,^,rs^y) 



5{t - Ty)5{xt 



(42) 



ys) 



in which rs.x{y) is as defined in Eq. (35). The same equa- 
tion appears in CC. We shall also suppress the subscript 
S unless necessary, as it is normally already specified in 
the notation for Go,b. 

The transverse part of the integral in Cg{S; n) can be 
factored out as 



jd-l 



xs±S'^ ^{xs± - ys±) 



Recalling our discussion about Eq. (11), this is a coarse- 
grained ^-function. At short distances, instead of a di- 
vergence, we should use 



s:d~i 



{xs± - ys±) 



iA/{2n)) 



(43) 



Therefore, the transverse direction integral becomes 
(A/(27r))'^-i / d"^-^- - 



xsi_ = iA/{27r)r-' f dSx ■ ns. 

JdA 

Identifying A'^^^ng as the surface element dSk, for 
a given patch the integration can be rewritten as 
{2tt)~'''~^^ <fg^\dSx ■ dSk\- This leaves us with only an 

integral over (g[)"^''(5; r^r, r'a;) - nG'^Ql{S;rx,rx)). 

The solution for the (1 -I- l)d Green's function on the 
n— sheeted replica manifold is given in CC: 



^ 00 

Goii^'^'^x^ry) = — dkCk/n{Ox - Oy) 



k=0 



(44) 



'X9k/7i{'rx,ry), 



where = l,(ij. = 2 for /c > 0, Ci,{9) = cos{y9), 
g,{r,r') = e{r - r')I,{r')K,{r) + e{r' - r)I,{r)K,{r'), 



11 



and /y(r) and K^{r) are the modified Bessel functions of 
the first and second kind respectively, r and 6 are again 
the polar coordinates of the ng — f plane, and we have 
suppress the index S' of r as only one patch direction is 
involved. 

The integral over G^^^ is 

■' k 

(45) 



So we have 



CGiS;n)-nCGiS;l) 



24n 



(48) 



Combining the above results into Eq. (40) and converting 
the sum over S into an integral around the Fermi surface, 
we obtain Eq. (19) for this geometry. 



The integral is divergent since the integrand 
rxgk/n{ra:,rx)\r^^oo= 1/4, a consequence of the 
fact that we are calculating the partition function of an 
infinite system. But this divergence should be canceled 
in CaiS^n) — nCaiS; 1). To regularize the divergence, 
we use the Euler-MacLaurin (E-M) summation formula 
following CC, and sum over k first: 



2 ^ poo 1 

-J2dkf{k)= f{k)dk--f{o) 

k=0 •'^ 



E 

i=2 



B 



(46) 



2i_ r(2j-l) 



(2j) 



(0), 



where i?2n are the Bernoulli numbers, /*^^^ ^^0) = 
i9^-'~^/(fc)|fc=o- Note that the first term, the integral over 
k, is always canceled by rescaling k/n ^ k in gk/n ■ For 
the remaining terms, which contain derivatives with re- 
spect to fc, we may add a constant, in this case —1/4, un- 
der the derivative, which allows us to pull the derivative 
outside the integral. The integrand now is well-behaved 
at infinity. To be more precise, according to Eq. (46), 
we need to compute 



D. Differential Equations of the Green's Functions 
and an Iterative Solution 



In this part, we derive the differential equations of the 
Green's functions for the quadratic boson theory with 
inter-patch coupling, and provide an iterative solution. 

Including the Fermi liquid interaction V{S,T;x ~ y) — 
UsT, the Hamiltonian becomes 



Y,9s,Tds<l^{S;x)dT<j^{T-x)\ 



(49) 



S,T 



drxrxdlgk/n{rx,rx) 



fc-j-O 



dl J drx(rxgk/n{rx,rx) - 1/4) 





(47) 



where gs.T — ^'^'J'*^ is order 1/A^. This Hamiltonian can 
be written in terms of the non-chiral fields as 



H ^ ^ d'^[Y. ii9-sV(S; x)f + (dxsXiS; x))^) + ^ gs,T{d.sX{S; x)dx^x{T; x) + 9^,^(5; x)dxMT: x))) , 

S S.T 

(50) 

I 



where we have made use of the fact that gs.T — g-s.-T, 
which is required by time-reversal symmetry. Here 
the summation over S is restricted to a semicircle. 
This Hamiltonian contains generalized type kinetic terms 
(inter-patch coupling due to interaction) which are not 
diagonal. To obtain the corresponding Lagrangian, one 



needs to invoke the general Legendre transformation[4.3], 
and obtains the following Lagrangian densities, respec- 
tively, in terms of ip or x- 
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s 



(h2iS, T)dMS; x)dMT; x) - A (5, T)ds^iS; x)dTv{T- x) 

S,T 

^ (h,{S, T)dtx{S; x)dtx{T; x) ~ h{S, T)dsx{S; x)dTx{T; x) 



S.T 



(51) 



where 

fi{S, T) = gs,T + 9-s-T — gs-T — 9-s,T, 
h{S, T) = gs,T + 9-S.-T + g-s,T + gs-T, 

and /ii(2)(S, T) is defined through 

{/+[A(2)(5,T)]ri = /+[/ii(2)(5,T)]. 



(52) 



Here, / is the identity matrix, and [f{h)i{S,T)] is the 
matrix formed by f{h)i{S,T), i = 1,2. Applying this 
result and making use equations of motion obtained from 
the Hamiltonian, we obtain the Lagrangians or C^. 
Here we arbitrarily choose to work with dp. Then, by 
making use of the E-L equation of motion, applying the 
same rescahng Eq. (37), and letting t = it, we obtain the 
differential equations that the interacting Green's func- 
tion G(") = {ip{S]x)(p{T]x')Ys satisfies: 

-{dl + dl-l)G^^\S,T-x,x') 

+ Y.^h2{i, T)dl + /i(;, T)a,aT)G(") (i, t- x, x') (53) 

I 

= CSs^rSir - t')5{xs - Xs))S{xs± - xsj_). 

Here the Jacobian due to change of variables is the same 
as in the free fermion case. The entanglement entropy 
is still given by Eq. (40), but replacing Gg"' with G'") 
in CaiSin). In the following, we omit the replica index 
n in the Green's function unless different values of n are 
involved in a single equation. 

As is well known, differential equations such as the 
above can be converted to an integral form[ i-t] relating 
the full Green's function to the noninteracting one. This 
leads to an iterative (perturbative) definition of the for- 
mer in terms of the latter. In the present case, this inte- 
gral equation reads 



G{S,T;x,y) 
= GoiS,T;x,y) + 

X (5^(/i2(/,T)a2 



d^zGo{S,S;x,z) 
h{l,T)didT)G{l,T;z,y) 



(54) 



GoiS,T;x,y)+SGiS,T;x,y). 



Given this equation, we can now compute the Green's 
function and thus the entanglement entropy perturba- 
tively in powers of U. 



E. Entanglement Entropy from the Iterative 
Solution 



In Eq. (54), the Go term is the same as that of the 
free fermions, thus yields the same contribution to en- 
tanglement entropy. To study how the correction term 
SG{S, T; X, y) affects the entanglement entropy, we need 
to study 



(fx5G{S,S]x,x)^ j d^x5^^'^G{S,S;x,x), 

M=l •' 

(55) 

where 5^^^'G denotes the Mth order correction. There 
are two distinctive types of terms in the perturbative 
expansion of 5G. In general, at order M , we have in total 
3(A/+1) integrals. Let us examine one of the many terms 
contributing to the M-th order correction, to be summed 
over patch indices: 



j d^xS^'^'^G{S,S;x,x) 

. M-l 

^ I d^xll {d^z,)Go{S,S;x,zo) (55) 

i=0 

X 9^^Go(Zo, ^0; ^0, ^i) ■ • • X dl.Goili,h', -Zi, ^li+i) 
X ■■■ X d'^^^_^Go{S,S;zM-i,x). 



Here we only include the r— derivatives. In general we 
would also have spatial (n;.) derivative terms, as well as 
terms with mixed derivatives. But f and ris directions 
are equivalent. Using rotational symmetry, and the fact 
that the two different derivatives in each term are with 
respect to independent variable that are each integrated 
over, one can see that all terms are identical except for 
S'-dependent pre-factors. The two categories of terms are 
defined by the set {h}: 1) h = S ^ i, i.e. with intra-patch 
coupling only; and 2) 3 li S containing inter-patch 
coupling. We shall label the two categories as 



6^^'^^G{S,S;x,x) = 

^LtL^i^' ^) + '^intor 



G{S,S;x,x) 



(57) 
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1. Intra-patch coupling and comparison with ID 

Setting li = S for all «'s in Eq. (56), first we consider 
the transverse direction 

Go{S, S; z,,z,+^) ^ 5{7}^l - z^f^,). 

We can immediately integrate out the transverse compo- 
nent of all Zi's and obtain 

/ dxs±S{xs± - zs±fl)S{zs±,M-i - ys±)\v^:B 

frr TT 

X / [[dzs±,i [[S{zs±.i - zs±,i+i) 

i 



In the last line we use again the fact that the transverse 
(5-function is a coarse-grained one (Eq.(43)). 

The rest of S^^^^^G is obtained by substituting 
Go{S, S; Zi, Zi+i) with the (1 + l)d Green's function 
Go^biS; zs^i, zs^i+i). Although a direct computation is 
possible, we first give a general argument that for any M 
the contribution to entanglement entropy from 6^^^j.^G 
vanishes. We do so by making a comparison with the ID 
case where a rigorous solution is available. 

For the ID Luttinger liquid with only forward scatter- 
ing, the entanglement entropy can be calculated directly 
via bosonization and the result remains at 1 /3 log L in 
the presence of interactions. The calculation is possible 
because, in our language, there are only two patches, so 
the transformation which diagonalizes the Hamiltonian is 
not plagued by the nonlocality issue we encounter in the 
2D theory. However, we can also treat the ID case with 
our perturbative approach. The resulting series of inte- 
grals turns out to be identical to the one obtained from 
the intra-patch contributions in the higher dimensional 
case except for the transverse (5— function. Therefore, we 
argue that at all orders, the intra-patch coupling terms 
have vanishing contribution to the entanglement entropy. 
We shall demonstrate such behavior explicitly up to sec- 
ond order in U later on. 



2. Scaling analysis of inter-patch coupling 

For terms with inter-patch coupling, wc find that they 
are of order 0{1/L) comparing to the leading term ac- 
cording a scaling argument. The crucial observation here 
is that, as long as 3 ^ S, we do not encounter the factor 
S°-\0) = L^-i(A)^-i, Eq. (43) because for ? 7^ 5 



J(.S)-.f))^(z«-.S) = ^ 



where Og (Oi) is the angle between ns {"hi) and the x- 
axis in the x — y plane. Therefore, when we integrate out 
the (M + 1) transverse (5-functions, the factor 6^^^{0) — 
L^^^(A)^^^ would be suppressed by even a single k ^ 
S. 

To examine the remaining integral, we can ignore the 
angular part as it cannot affect the scaling behavior. The 
asymptotic expansion of K^(r) and /^(r) for real r at 
large value is [45] 



lAr) = 



/27rr 



1+ E 

n=l 



(2r)" 



where {iy,n) = Jr(i/2T.^"«) ■ ^y using the 
above asymptotic expansion of Bessel functions, 
the leading term for d'^^Go{li,l-Zi, Zi+i) behaves as 

~ 9iz, - z,+i)e-(^'-^'+i)/(z, - z,+i) + 9{z,+i - 
Zi)e~(^'+i~^'^ /{zi^i — Zi). All of these terms peak around 
Zi+i = Zi and are otherwise exponentially suppressed. 
We may therefore again estimate this integral by letting 
X = Zq = Zi ~ ■ ■ ■ — zm-1 and removing (M -I- 1) of the 
integrals. The remaining integrals yield, at the leading 

However, at the 



order, /d^^+^zl/; 



M 



J dzz 



M ^-M 



leading order, there is no v dependence. According to 
the formalism in Sec. IV C, such terms have no contri- 
bution to the entanglement entropy. Therefore, the term 
that contributes to the entanglement entropy is the next 
order which behaves as J dz^ and is of order ©(logL), 
leading only to a correction ~ O(logL) x logL to the 
entanglement entropy. 

Next, we shall demonstrate in detail our above analy- 
sis, for both inter-patch and intra-patch coupling terms 
by explicit calculation up to the second order. 



3. First Order Correction 



The first order term correction to J dxG{S, S; x, x) is 



6^^'>CG{S;n)^ J d^xd^ zGo{S, S; x, z) 
X {h2{S,S)dl+fi{S,S)djJGo{S,S;z,x). 



(60) 



As we have pointed, it is sufficient to calculate either 
piece of the two terms due to the equivalence of the imag- 
inary time direction and the real space direction. The 
other piece should be just the same except for the coeffi- 
cient. Here we choose to compute the first term. 

The transverse degrees of freedom provide an overall 
factor counting the total degrees of freedom as discussed 
in the general case. Then we can also integrate out the 
angular degrees of freedom in the — f plane, both 9x 
and 9z as defined in Eq. (44), after which one obtains 



smit/; 



(59) 



S^'^CaiS; n) ^ V ^S^'^G^/n <f {dS, ■ dS. 



(61) 
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where 



(62) 



X id' 



:)9k/n{rx,r^). 



The two k summation is reduced to one due to orthogo- 
nahty of the angular function Cfe/„(0). By employing the 
E-M formula and properties of the Bessel functions, we 
show in Appendix A that sum over fc-values in Eq. (61) 
can be converted into an integral, which cancels in Eq. 
(40) for the same scaling reasons discussed above, follow- 
ing Eq. (46). Therefore, we find that the contribution of 
Eq. (60) to the entanglement entropy vanishes. 



4- Second Order Correction 
The second order correction is 
(5(2)CG(S';n) = j d^xd^zd^ ziGo{S , S;x, z) 
xY,{h2{l,S)dl+fi{l,S)d,,d,,')G„{l,l;z,zi) (63) 
X (h2{S,S)dl^^+h{S,S)dUGo{S,S:z,,x). 
• for Z = S: 

J d^xd^zd^ziGo{S,S;x,z) 

X {h2{S,S)dl+h{S,S)dl)Go{S,S;z,z^) (^4) 
X {h^iS, S)dl^ + his, S)dlJGoiS, S; z,,x). 

According to our general discussion, we only need to con- 
sider the following piece: 



After a lengthy but similar calculation as for the first 
order (see Appendix B), we find, using the E-M formula: 



dr^dr^dvir^r^ri 
dkpk/n{rx,rz,ri) 



(67) 
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n 

16fc' 



where Pk/ni'i'xjf'zTri) is the product of g^./^ dependent 
terms in Eq. (66). The usual scaling argument for the 
integral shows that the entire expression is proportional 
to n, and thus cancels the second {n = 1) term in Eq. 
(40): 



S{S) ^ - 



d_ 
dn 



d'^x{Gn-nGi)) 



Therefore, at the second order level for the / = S piece we 
still have no correction to the scaling law of entanglement 
entropy. 

. for Z ^ S: 

The integrand we need to consider is 
GoiS,S;x,z) 

X (^h2{l,S)dl^+h{l,S)dz,dz,)Go{l,l;z,zi) (68) 
X {h2{S,S)dl^^ +h{S,S)dlJGo{S,S;zi,x). 

The first thing to notice in Eq. (68) is that we have 
derivatives along directions different from the patch nor- 
mal direction ng acting on the non-interacting Green's 
function. We expand this term as 

dzidzsGoil,l;z,zi) Siz^f' - z[^l)dzidzsGo,bil; z, zi) 



+ dzJizf' - z['l)dz,G„.b{l;z,zi). 



(69) 



For the first term, we can decompose the derivative 



/ d^xd^zd^ziGa{S,S\x,z)dlfio{S,S;z,Zi) into terms that act along fii and along its trans- 

J verse direction, resnectivelv- The riou-intera.ctirie' Oreeu'f 



X dl^ Go{S,S;zi,x) 



(65) 



where 

<5^^-'Gfc/„ = J dr^dr^dri {rx,rz) 

X (dl-ky{rznr)g,/^{rz,r,) ^^^^ 
X {d^^ - k^ /{rinf)gk/n{ri,r^). 

In the above, we have proceeded as in the first order cal- 
culation, integrating out the angular part first to obtain 
the expression for (5'^^'Gfc/„. 



verse direction, respectively. The non-interacting Green's 
function only depends on the transverse coordinates 
via Go{l,l]x,y) ^ d{x''^ — yf'), which indicates that 
those derivative terms vanish. Thus it is ^ S{z^f^ — 

Zi±)dz,Go_bil', z, Zi). For the second term, we integrate 
by parts with respect to zg, which leads to (including 
now the first Go factor, which depends on z) 

-S{xf'>-z[l^)6{zf^-z[l)dzsGoAS; X, z)dz,Go,b{l; z, z,). 



Therefore, the overall integrand is proportional to 

6{x^f'' - - x^f'>)Siz^^^ - zfl). Note that the 

(s) 

Cj_ dependence only appears in these (5-functions, we 



can integrate it out, leaving only S(z^^^ — z[^)S{z 



(0 



41) - ^(^1 
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Secondly, it is sufficient to focus on the following terms 
in the integrand 

(GQ^b{S;x, z)dl^GQ^b{l; z, zi) + d,^sGo,biS;x, z) 

X d^fiQ^i,[l;z,Zi)]dl Go,b{S;zi,x), 



to ease the presentation. For other combinations, the rest 
of this section is equally applicable with minor modifica- 
tions that only leads to different coefficients and do not 
affects the scaling analysis. We first perform the intra- 
patch integration 

J d^xGo,biS; X, z)Go,biS; z,,x) = H{S; z, z,), (70) 

where 

H{S;z,z,)=Y, ^Ck/n {d. , e,, ) {e{r, - r, ) 

A; 

So for a given S the contribution to entanglement en- 
tropy due to coupling with patch I can be written as 



1 



J d^zSzi5^{zi-z)(dl^^H{S;z,zi) 



X d^Ga^b{l]Zi,z) + d^sd^^^H{S; z, zi)dz,Go,b{l; zi, z) 
1 



sin(0i - Os) 

|2 



J d?zdT^ (dl^^ H{S- z, z^)dlGo,b{l\ z^,z) 



^zs^r,, H{S\z, Zi)d^^Go^b{l-i zi,z) 



(71) 



where z^x,z^y indicate the two spatial components of z. 

As we argued in previous section, for extracting the 
order of magnitude of the result it is sufficient to set 
Ti = r in the final line of Eq. (71) and remove the 
integral over ti. We also note that the derivatives do 
not alter the leading power of r, owing to the presence 
of the exponential function. Therefore, it is sufhcient to 
examine 



J £z{H{S;z,Zi)Go,b{l\z^,z)) 
At the lowest order in 1/r, we have 

H{S-r,r) ^ rI,{r)K,+ {r) - rK,{r)h^{r) 



(72) 



(73) 



=(1 



(1 



(j.-l,l) 



)(1 



2r 



2r 



)(1- 



(M) 
2r 



+ ...) 



(74) 



Since the r derivative does not alter the leading powers, 
we extract the leading term to be 



[H[S-z,zi)dlGo,b{S;zi,z)) 



1 



(75) 



For a triple integral over l/z^, one would get a linear 
divergence, i.e. the result would be '-^ L. This is indeed 
the case as we have already seen in previous calculation. 
However, at the lowest order, everything is independent 
on V ~ k/n. Actually what finally appears in the the en- 
tanglement entropy are the fc— derivatives of these terms 
appearing in the E-L summation formula. This means 
the leading term has vanishing contribution to the en- 
tanglement entropy. The first term contributing to en- 
tanglement entropy is then ~ J <Pz\ the upper limit of 
which is order OilogL) and only leads a correction up 
to ~ ©(logL) X logL to the free fermion entanglement 
entropy. 



V. SUMMARY AND CONCLUDING REMARKS 

In this paper, we developed an intuitive understand- 
ing of the logarithmic correction to the area law for the 
entanglement entropy of free fcrmions in one and higher 
dimensions on equal footing - the criticality associated 
with the Fermi surface (or points). Then we used the 
tool of high dimensional bosonization to compute the en- 
tanglement entropy, and generalized this procedure to 
include Fermi liquid interactions. In the presence of such 
interactions we calculated the entanglement entropy for 
a special geometry perturbatively in powers of the in- 
teraction strength up to the second order, and find no 
correction to the leading scaling behavior. We also point 
out that the situation is the same at higher orders. Our 
results thus strongly suggest that the leading scaling be- 
havior of the block entanglement entropy of a Fermi liq- 
uid is the same as that of a free Fermi gas with the same 
Fermi surface, not only for the special block geometry 
studied in this paper, but for arbitrary geometries. Ex- 
plicit demonstration of the latter is an obvious direction 
for future work. 

In the special geometry in which we performed explicit 
calculations using the replica trick, a mass- like term is in- 
troduced to regularize the theory at long distance, as is 
done in closely related contexts[. ., ->-]. For a Fermi liq- 
uid (which is quantum-critical) the corresponding length 
scale ^ ~ vp/m must be identified with the block size 
L, and is thus not an independent length scale. On the 
other hand, such a mass-like term can also describe a 
superconducting gap due to pairing. In particular, for 
a weak-coupling superconductor, ^, the superconducting 
coherence length, is much longer than all microscopic 
length scales, but finite nevertheless. In this case it is 
independent of L, and the interplay between the two is 
interesting. For L < ^, the Fermi liquid result (1) still 
holds. But for L > ^, the logarithmic factor in the en- 
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tanglement entropy saturates at log^, and we expect: 

SiPA) = T979^ loge X / / \dS^. -dSkl (76) 
12(27r)^ JoAJdr 

which agrees with the conjecture made in Ref.['ii^]. 

More generaUy, Fermi hquids are (perhaps the best un- 
derstood) examples of quantum critical phases (or points) 
in high dimensions. Unlike in ID where conformal sym- 
metry powerfully constrains the behavior of entangle- 
ment entropy, our understanding of entanglement prop- 
erties of such high-dimensional quantum critical phases 
or points (many of them have Fermi surfaces but are not 
Fermi liquids) is very limited. Our work can be viewed 
as a step in that general direction. Furthermore, the for- 
malism developed in this work has potential applicabil- 
ity to systems with composite or emergent fermions with 
Fermi surfaces as well, or more generally, non-Fermi liq- 
uid phases with Fermi surfaces. The system studied in 
Ref. [^"'2], where there is an emergent spinon Fermi sur- 
face, is a potential example. 
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Appendix A: Calculation of (5^^'Gfc/„ 

Throughout the Appendix, we shall denote the modi- 
fied Bessel functions K^{ri), I^{ri) as Ki, li for simplicity 
with V = k/n. We also have K{I)y±i{ri) which shall be 
shortened as K[I)i^±. 

i5^^'Gfc/„ = / dr^dr^r^r^gk/n{r^,r^) 

(Al) 



)9k/n{rz,rx). 



Expanding (9^^ - :;^)gk/nirz,rx), and noting the iden- 
tities I'K - K'l =\/x, X" - {v^/x^)X = X- {l/x)X' 
where X = K,I, we get 

= - ^ii:^— ^ + eir, - r,)(l - -drJIzK^ (A2) 

+ 0ir,~rx){l- -dr.)hKz. 

Then integrating over r^ first, and making use of the 
following formula 

/ dxxXl{x) = ^^x'iXlix) - X,_i(x)X,+i(x)), (A3) 



where Xyix) can be the first or second kind of modi- 
fied Bessel function, ly or K^, and the identities l^K + 
IK_ = I^K + IK^ = 1/x in addition to those given 
above, S'^^^G^/n is reduced to 

S^^^Gk/n^ J drr(^-IK + r^{I^ ~I+I^)K^y (A4) 

We apply the same strategy as in Sec. IV C, making 
use of the E-M formula to do the sum over the fc-index 
in Eq. (61). This converts the sum into a divergent 
integral over k which cancels in Eq. (40) as before, and a 
sum over terms of the form d^S^^"^ G k / n\ /.^q that turn out 
to vanish, as we will now show. Again, we can include 
proper constants under the derivative into the integrand. 
These derivatives then act on well defined integrals. The 
first term has been discussed in Sec. IV C: 



(A5) 



drdl{-rIk/nKk/n) 
= dildri~rIK+\)=di{^). 
The second term can be shown to be 
J drdlir^I^ - I+I^)K^) 
= dij dr{r\l^-I+I^)K^-\) 



(A6) 



^^^( 16 2n 



Summing the two terms together, we find 



G 



k/ n 



(A7) 



for all J > 0. 



Appendix B: Calculation of S^-^^'Gk/n 

Let us first compute the integral: 

dr2r2gi,{r,r2){dl^gy{r2,ri) - 
X gk/nir2,ri)j 

= 9{r — ri)h{r, ri) -\- Oiji — r)h{ri, r) 



(Bl) 



with 



/i(r,ri) = -/i/v +^(/i(r,ri) + i^/iln-), 

/i(r,ri) = KK^ill - h.,+h.-) - I'D - Ih 
X {r'^iK'^ - K+K_) - K^) + Kh{F{r) - F{ri) 
-IK + hK,), 

F{r) = 2 / drrIK = r^IK + r^I+K_. 
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The IKi term in h{r, ri) results from (5-functions (deriva- 
tives of (5-functions) coming from the derivative applied 
on the step function (9{r)). The remaining part comes 
from terms involving a product of two ^-functions. Here 
one needs to distinguish between r > ri and r < ri, 
which gives rise to the terms in 0{r — ri) and 9(ri — r), 
respectively. The remaining integral 



drdrirri{9{r — ri)h{r, ri) + 9{ri — r)h{ri, r)) 



X I - ^(^^^ + 9in - r){l - -dr, )KJ (B2) 
r ri 



+ 9{r-n){l- -dr,)Kh) 



ri 

can be carried out by applying the identities of Bessel 
functions / and K used in Appx. A. In applying the E-M 
formula to the sum over k in (65), we again arrive at a di- 
vergent k-integral that can be rescaled and subsequently 
canceled (see (67) and below), and a sum over deriva- 
tive terms that are well-behaved. In the latter terms, 
we always add proper constants under the derivatives to 
regularize the integrand at infinity, as before. We divide 
(B2) into two terms. The first is the one containing the 
(5-function. After integrating out ri , this term becomes 

dl{- / drrh{r,r)) ^ dl{~ / dr{rh{r,r) + ^)) 



(B3) 



= dl{~J dr{r{-IK + "^{-K^I+I^ 



The second term is expanded to 



drdrirri (j){r — ri)h{r, ri)K(Ii — I[/ri) 



(B4) 



Due to the complexity of h{r,ri), we examine each of 
the three terms of h{r,ri) separately. The first term is 



simple. Applying those identities of K's and /'s and in- 
cluding the proper constant, we get 



5^ J drdri{rri{9{r ~ ri){-hK)K{Ii - I[/ri) 
+ 9{r^ - r){-IK,)I{K, - K[/r,)) - J)) (B5) 



The second term of h{r,ri) contributes 

di[^ J drdnrn9{r-n)fl{r,n)K{h-^^I[) 

+ drdr,rn9in - r)/l(ri, r)(iCi - ^^K[)I 

(B6) 

By interchanging the dummy variables r and ri , employ- 
ing the properties of the modified Bessel functions with 
care, and including the regularization constant, we arrive 
at 



dl{l J drdri{rri9{r-ri)fi{r,n)i2Kh^KI[/n 



(B7) 



The last part of h{r, ri) contributes as 

J drdri^^9{r - ri)Kh\n^{K' h/r 



(B8) 



KI'Jt,) = - 



Summing all the above terms together we get 

a^,(5(^)G,/„)-a^(-^). (B9) 



[1] M. A. Nielsen and I. Chuang, Quantum Computation 
and Quantum Information, Cambridge University Press, 
Cambridge (2000). 

[2] L. Amico and R. Fazio, Entanglement and magnetic or- 
der, J. Phys. A: Math. Theor. 42, 504001 (2009). 

[3] J. I. Latorre and A. Riera, A short review on entangle- 
ment in quantum spin systems, J. Phys. A: Math. Theor., 
42: 504002 (2009). 

[4] J. Eisert, Colloquium: Area laws for the entanglement 
entropy. Rev. Mod. Phys., 82, 277 (2010). 

[5] S. R. White, Density matrix formulation for quan- 



tum renormalization groups, Phys. Rev. Lett., 69, 2863 
(1992). 

[6] F. Verstraete, D. Porras, and J. Cirac, Density Matrix 
Renormalization Croup and Periodic Boundary Condi- 
tions: A Quantum Information Perspective, Phys. Rev. 
Lett., 93, 227205 (2004). 

[7] D. Perez-Garcia, F. Verstraete, M. M. Wolf, and J. I. 
Cirac, Matrix Product State Representations, Quantum 
Inf. Comput. 7, 401 (2007). 

[8] M. Levin and X.-G. Wen, Detecting Topological Order 
in a Cround State Wave Function, Phys. Rev. Lett., 96, 



18 



110405 (2006). 

A. Kitaev and J. Preskill, Topological Entanglement En- 
tropy, Phys. Rev. Lett., 96, 110404 (2006). 
H. Li and F. D. M. Haldane, Entanglement Spectrum as 
a Generalization of Entanglement Entropy: Identification 
of Topological Order in Non-Ahelian Fractional Quantum 
Hall Effect States, Phys. Rev. Lett. 101, 010504 (2008). 
X. Chen, Z.-C. Gu, and X.-G. Wen, Classification of 
gapped symmetric phases in one- dimensional spin sys- 
tems, Phys. Rev. B 83, 035107 (2011). 
T. Osborne and M. Nielsen, Entanglement in a sim- 
ple quantum phase transition, Phys. Rev. A, 66, 032110 
(2002). 

G. Vidal, J. I. Latorre, E. Rico, and A. Kitaev, Entangle- 
ment tn Quantum Critical Phenomena, Phys. Rev. Lett., 
90, 227902 (2003). 

For a review, see J. Eisert, M. Cramer and M. B. Ple- 
nio. Colloquium: Area laws for the entanglement entropy. 
Rev. Mod. Phys. 82, 277 (2010). 

C. Holzhey, F. Larsen, and F. Wilczek, Geometric and 
Renormalized Entropy in Conformal Field Theory, NucL 
Phys. B, 424 (1994) 443-467. 

D. Gioev and I. Khch, Entanglement Entropy of 
Fermions in Any Dimension and the Widom Conjecture, 
Phys. Rev. Lett., 96, 100503 (2006). 
M. Wolf, Violation of the Entropic Area Law for 
Fermions, Phys. Rev. Lett., 96, 10404 (2006). 
M. Cramer, J. Eisert, and M. B. Plenio, Statistics De- 
pendence of the Entanglement Entropy, Phys. Rev. Lett., 
98, 220603 (2007). 

W.F. Li, L.T. Ding, R. Yu, T. Roscilde, and S. Haas, 
Scaling Behavior of Entanglement in Two- and Three- 
dimensional Free-fermion Systems, Phys. Rev. B, 74, 
073103 (2006). 

T. Barthel, M.-C. Chung, and U. SchoUwock, Entangle- 
ment Scaling in Critical Two-dimensional Fermionic and 
Bosomc Systems, Phys. Rev. A, 74, 022329 (2006). 
G. C. Levine and D. J. Miller, Zero Dimensional Area 
Law m a Gapless Fermionic System, Phys. Rev. B, 77, 
205119 (2008). 

Y. Zhang, T. Grover, and A. Vishwanath, Entanglement 
Entropy of Critical Spin Liquids, Phys. Rev. Lett., 107, 
067202 (2011). 

R. Helling, H. Leschke, and W. Spitzer, A Special Case of 
a Conjecture by Widom with Implications to Fermionic 
Entanglement Entropy, Int. Math. Res. Notices (2011) 
2011: 1451. 

A preliminary discussion of this idea, as well as its con- 
nection with multi-dimensional bosonization was pre- 
sented earlier by one of us[_ >]. A similar idea has been 
independently discussed by B. Swingle[26], based on a 
gapless mode-counting argument which is equivalent to 
our toy model. In this work, we construct the local de- 
grees of freedom explicitly via high-dimensional bosoniza- 
tion, and use the constructed bosonic theory to compute 
the entanglement entropy directly. We also generalize the 
calculation to include Fermi liquid interaction in a spe- 
cial geometry, which is non-trvial and shows subleading 
contribution ~ 0(1) x InL that cannot be accounted for 
by mode-counting. 
[25] W. X. Ding, Entanglement Entropy in Many-Fermion 
System, arXiv:0908.2418. 



[26 



[27 



[28 



[29 

[30 

[31 
[32 

[33 
[34 

[35 

[36 

[37 

[38 

[39 

[40 
[41 

[42; 

[43 
[44 

[45 



B. Swingle, Entanglement Entropy and the Fermi Sur- 
face, Phys. Rev. Lett. 105, 050502 (2010). 
B. Swingle, Conformal Field Theory on the Fermi Sur- 
face, arXiv:1002.4635vl, and Renyi entropy, mutual in- 
formation, and fluctuation properties of Fermi liquids, 
arXiv:1007.4825vl. 

F. D. M. Haldane, Luttmger's Theorem and Bosonization 
of the Fermi Surface, in Proceedings of the International 
School of Physics "Enrico Fermi", Course CXXI "Per- 
spectives in Many-Particle Physics", R. A. Broglia and 
J. R. Schneffer eds., North-Holland, Amsterdam (1994), 
pp 5-29 (1993); arXiv:cond-mat/0505529. 
A. Houghton and J. B. Marston, Bosonization and 
fermion liquids in dimensions greater than one, Phys. 
Rev. B, 48, 7790 (1993). 

A. Castro Neto and E. Fradkin, Bosonization of the low 
energy excitations of Fermi liquids, Phys. Rev. Lett., 72, 
1393 (1994a). 

A. Castro Neto and E. Fradkin, Bosonization of Fermi 
liquids, Phys. Rev. B, 49, 10877 (1994). 

A. Castro Neto and E. Fradkin, Exact solution of the 
Landau fixed point via bosonization, Phys. Rev. B, 51, 
4084 (1995). 

H. Widom, On a class of integral operators with discon- 
tinuous symbol, Oper. Th.: Adv. Appl, 4, 477 (1982). 
There is a recently announced proof of Widom's conjec- 
ture by A. V. Sobolev, Quasi- Classical Asymptotics for 
Pseudodifferential Operators with Discontinuous Sym- 
bols: Widom's Conjecture, Funkts. Anal. Prilozh., 44:4 
(2010), 86-90. 

V. E. Korepin, Universality of Entropy Scaling in One 
Dimensional Gapless Models, Phys. Rev. Lett., 92, 
096402 (2004). 

B. -Q. Jin and V. E. Korepin, Quantum Spin Chain, 
Toeplitz Determinants and the Fisher- Hartwig Conjec- 
ture, J. Stat. Phys., 116, 79 (2004). 

P. Calabrese, J. Cardy, and E. Tonni, Entanglement en- 
tropy of two disjoint intervals in conformal field theory, 
J. Stat. Mech.: Theory Exp., 0911, PllOOl (2009). 
V. Alba, L. Tagliacozzo, and P. Calabrese, Entanglement 
entropy of two disjoint intervals in c=l theories, J. Stat. 
Mech.: Theory Exp., 2011, P06012 (2011). 
S. Furukawa, V. Pasquier, and J. Shiraishi, Mutual In- 
formation and Boson Radius in a c=l Critical System in 
One Dimension, Phys. Rev. Lett., 102, 170602 (2009). 
R. Shankar, Renormalization- group approach to interact- 
ing fermions. Rev. Mod. Phys., 66, 129 (1994). 
P. Calabrese and J. Cardy, Entanglement entropy and 
quantum field theory, J. Stat. Mech.: Theory Exp., 2004, 
P06002 (2004). 

M. Hertzberg and F. Wilczek, Some Calculable Contri- 
butions to Entanglement Entropy, Phys. Rev. Lett. 106, 
050404 (2011). 

H. Goldstein and C. P. Poole, Classical Mechanics, Ad- 
dison Wesley (2001). 

S. Doniach and E. H. Sondheimer, Green's Functions 
for Solid State Physicists, World Scientific Publishing 
Company (1998). 

Z. X. Wang and D. Guo, An Introduction to Special 
Function, Peking University Press (2000). 



